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ABSTRACT 

Physical solutions to convex scalar conservation laws satisfy a one-sided 
Lipschitz condition (OSLC) that enforces both the entropy condition and their 
variation boundedness. Consistency with this condition is therefore desirable 
for a numerical scheme and was proved for both the Godunov and the Lax- 
Friedrichs scheme — also, in a weakened version, for the Roe scheme, all of 
them being only first order accurate. A new, fully second order scheme is 
introduced here, which is consistent with the OSLC. The modified equation is 
considered and shows interesting features. Another second order scheme is 
then considered and numerical results are discussed. 
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1. The OSLC for Convex Conservation Laws. For a given (say locallv 

integrable) function v(x), it makes sense to consider the one-sided Lipschitz semi- 
norm: 

(1.1) p(u) = supf^^ — where a + = max(0,a)). 

xjty ' x ~ y 

If v is a decreasing function, then p(v) = 0. Conversely, if v has increasing jumps, 
then p(u) = + 00 . 

It has been known for a long time [7], that the physical solutions to a scalar 
conservation law: 


(1.2) u t + f[u) x — 0 
with a (strictly) convex flux function: 

(1.3) /"(«)> o>0, 


satisfy: 


(1.4) 


„«,*) -y(t.y), J_ , for ( > 0 and X / s,). 
x — y at 


that is, in terms of the semi-norm: 
(1.5) 


p(u(f, •)) < — 7 for t > 0. 

at 


In particular, the semi-norm of the solution at time t > 0 is always finite, even 
when the semi-norm of the initial data is infinite. In other words, any increasing 
discontinuity introduced in the initial data is immediately spread out and becomes 
a rarefaction wave - and then only decreasing jumps, or shock waves, can appear. 
This is precisely the correct entropy condition for convex conservation laws. It is 
therefore automatically enforced by the OSLC (1.5). 


Moreover, the total variation boundedness of the solutions (a well known prop- 
erty of any scalar conservation law) is also enforced by the OSLC (1.5): Let us 
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consider, for simplicity, the initial value problem associated with (1.9), with peri- 
odic boundary conditions: 

(1.6) tt(<,0) =u(t,L)yt >0, 
where L > 0 is the period. Then the following estimate 

(1.7) TV{u(t, •)) < 2 Lp(u(t, •)) < 2 L(at)- 1 
holds, as a consequence of the following lemma: 

LEMMA 1. Let u(x) be a locally integrable function of period L, then 

(1.8) TV{v) < 2p(v)L. 


Proof. The proof is given when v is smooth and can be easily generalized by a 
standard density argument. We have: 

TV(v)= f \v'{x)\dx < f (|t/(x) - p(v)| + p(v)}dx 

Jo Jo 


By the definition of p(v), we have: t/(x) < p(v) and thus |u'(x)— p(t’)| — p(v) — v'(x). 
Therefore: TV (v) < / Q L [2 p(v) — u'(x)]dx = 2Lp(v) - v'{x)dx = 2 Lp{v) (since v 
is periodic). Q.E.D. 


So, at least for the periodic boundary problem, (in fact, also for the Cauchy 
problem), the OSLC enforces both entropy consistency and total variation bound- 
edness. 


As stated in (1.5), the OSLC is not entirely satisfying, since (a£) -1 blows up 
when t approaches 0. A more refined estimate can be quite easily deduced for 
p[u(t, •)) through the viscosity method: 
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PROPOSITION I. Any physical solution u(t,x) to (1.2-3) satisfies: 

(1.9) p(w(<, •)) < (p(«(s, -)) -1 + or(f ~ s)) -1 for t>s> 0 
and in particular, 

(1.10) p(u(t, •)) < p(u(v)), for t>s> 0. 

We do not- claim that this is a new result. 

Before proving this proposition, let us notice that, when p(u(0, •)) is finite (in 
particular when u(0, x) is smooth), then p{u(t , •)) is uniformly bounded in time. 

Proof (sketch) Let u e (t, x) be the (presumably unique and smooth) solution to: 

( 1 - 11 ) < + /«) = <*• 

A trivial computation leads to: 

(1.12) ^/>K(L ')) - )) 

that is (after integration) 

(1.13) p(u e (<, •)) < (p(u e {s, -)) -1 + (*{t — s)) _1 for t > s > 0. 

Then (1.10) is obtained when e -* 0. Q.E.D. 

II. Convergence of OSLC consistent numerical schemes. As usual, an 
explicit conservative numerical scheme is defined by: 

(2.1) « +1 -»")Ai+(/” +1 -/?_ 1 )A(=0 

' 2 2 

where A = (A t, Ax) denotes the time and space steps, u'f is an approximation 
to u(nAt,i£x) and f? +x is the numerical flux satisfying the usual requirements 
fi+x(ui+k, ..., tif—k+i ) is a Lipschitz function of 2k variables with 


fi+i(u,u,...,u) = /(«). 
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Let us introduce: 

j 

(2.2| P "=Pi=sup(^2^) + 

j 

DEFINITION. A numerical scheme is said to be OSLO (respectively weakly 
OSLO) consistent if: 

(2.3) < (p(u( 0, -)) _1 + cmAf) -1 
respectively: 

(2.4) Pa<p(«(0,-)) 

Remark: Conditions (2.3) and (2.4) clearly are discrete versions of (1.9) and 

( 1 . 10 ). 

We have the following result: 

PROPOSITION 2. For the periodic boundary initial value problem, the approx- 
imate solutions converge to the correct entropy solution , if either (i) the scheme is 
OSLC consistent, the initial condition has a finite number of increasing jumps, and 
the propagation speed of the scheme stays uniformly bounded; or (ii) the scheme is 
weakly OSLC and the semi-norm of the inital data p(u(0, •)) is finite. 

Basically the proof follows from the fact that the OSLC condition enforces both 
entropy consistency and total variation boundedness. Technical details are given in 
Appendix 1. 

Remark. The same result holds for the Cauchy problem when the initial value 
u(0, x) is constant for large values of |x|. 
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III. OSLC consistency of some first order accurate schemes. Properties 
(2.3) have been obtained in [4, 7, 12] for the Godunov and the Lax-Friedrichs 
scheme. It is unclear, to the best of our knowledge, whether or not this property 
holds for any monotone scheme (or more generally E-scheme, as defined in [8]) It is 
also clear that the Engquist-Osher scheme is OSLC consistent (as a matter of fact, 
the EO scheme is the discrete version of the transport-collapse method described 
in [2], which is clearly OSLC consistent). A more recent result involves the Roe 
scheme [3]: 

PROPOSITION 3. [3]. The Roe scheme is weakly OSLC consistent. When the 
initial condition satisfies p(u(0, •)) < -f-oo, in particular, when it is smooth, then 
the approximate solutions converge to the correct physical solution. 

This is a little surprising, since the Roe scheme is usually quoted to be entropy 
inconsistent! As a matter of fact, the Roe scheme does not violate the entropy 
condition, provided the semi-norm of the initial data, p(u( 0 , •)) is finite. Otherwise, 
in particular when the initial data has increasing jumps, the entropy condition can 
be violated, which explains the apparent contradiction. 

IV. MUSCL schemes and the OSLC. A priori, it is unclear whether the 
OSLC is compatible with second, or higher, order accuracy. For example, the E- 
condition, introduced in [8] to enforce both entropy consistency and total variation 
boundedness, for any kind of scalar conservation laws (not necessarily convex), is 
known to be only compatible with 1st order accuracy [8]. An appealing framework 
in which to discuss higher order accurate schemes, is the one introduced by Van 
Leer [13] with his MUSCL schemes, a generalization of Godunov’s ideas. In a rather 
abstract way, a MUSCL scheme can be described in terms of the exact evolution (or 
solution) operator (E t t > 0) associated with (1.2) and a projection operator P A * 
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that transforms any given function t>(z) into a cell wise polynomial approximation 
F Al r(x), in such a way that, in each cell, the mean value is preserved: 

P&xV{x) = Vi 



1 


(4.1) 

Ax 

■~2 

where 

_ Aii,. , 


(4.2) 

v(x)dx 


The numerical scheme is then defined by: 


(4.3) 


1 /’ x ' + i 

< = rj x ;«■<*>*• 

1 o 


where the function u n (x) is given by: 

(4.4) u n = (E± t P* x ) n u( 0,) 

The Godunov scheme corresponds to the simplest choice of the projection operator 

P Ax- 

Paxv{x) = Vi, for x,_A < x < x i+ 1 , where Vi is given by (4.2). 

This scheme uses piecewise constant approximations and is only first order 
accurate. Therefore Van Leer’s idea was to increase the order of accuracy by using 
polynomial approximations of higher order, when choosing Pax- 

In order to keep stability and avoid spurious oscillations, Pax is usually required 
not to increase the total variation: 


(4.6) 


TV(Paxv) <TV(v). 


Van Leer found a second order operator Pax, using piecewise linear approxi- 
mations and still compatible with (4.6). In each cell, 


(4.7) 


Paxv(x) = Vi + (x - Xi)si where z,- = iAx, 
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and the numerical slope s,- is given by 


(4-8) 

with 



Vj+i ~ Vi 
Ax 


Vi - Vj-l \ 
Ax /’ 


{ 1 if a, b > 0 

— 1 if a, b < 0 

0 otherwise 


Entropy consistency of the resulting scheme is a priori unclear and was proved 
in [9] (for the method of lines, i.e. the At — ► 0 limit case) and in [10] (provided 
further minor restrictions are added). This scheme is not OSLO consistent, mainly 
because the corresponding projection operator may increase the one-sided Lipschitz 
semi-norm p. Conversely we have: 


PROPOSITION 4. A MUSCL scheme designed in such a way that: 
(4.10) p{P& x v)<p{v) 

is automatically OSLO and therefore convergent. 


Proof. The proof is based on the fact that the exact evolution operator E t 
satisfies: 

(4.11) p{E&tv) < (p(v) -1 + <*A/) -1 

as a consequence of (1.9). 

Therefore, when (4.10) holds, we get: 

p(u n ) = p[E& t P&xU n ~ l ) by definition (4.4)) 

< {p{P& x u n ~ x )~ l + aAf) -1 (by (4.11)), 

< {pin"- 1 )- 1 + a At)- 1 (by (4.10)), 

< • • < (p(«(0, ))~ l + anAt) -1 (by induction). 
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Then, since, for each t, 


u? +l - «? _ _1_ /* Ax ti n (x, + i + a) -M n (x,_i +a) 


Ax Ax ; 0 

< p(u n ) (by definition of p), 

it follows: 


Ax 


ds by definition (4.3)} 


(4.12) PI = < P{u n ) < (P(«(0, •)) -r atn&t) 

which is exactly the discrete OSLO (2.3). This completes the proof. 


-l 


V. A fully second order OSLO consistent MUSCL scheme. Condition (4.10) is 
not satisfied by the classical MUSCL scheme corresponding to (4.7. 8, 9). Never- 
theless. 


PROPOSITION 5. The MUSCL scheme defined by (4-7,8) and: 

(5.1) o(a, b) = max(a, b) 

is fully second order accurate, satisfies (4-10), therefore is OSLO consistent and 
convergent. 

Proof. Let us first prove (4.10). We have to show, that for any x and y, x ^ y, 
we have 

(5.2) -A ; y . <y> < pW . 

x-y 

Because of definition (4.7), it is enough to check both 

(5.3) s, < p(v) and 

d e f 

(5-4) b i+ i = P Ax v{x i+ i +0) - P Ax v{x i+ i -0) < 0, 

hold for each i. [In other words, all discontinuities are decreasing and, in each cell, 
the slope is bounded from above by p(v)]. 
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Let us check (5.3). We have (by definitions (4.8) and (5.1)) 


(Vi 

Si = maxi — 


Vi + i - Vi 


Ax 


= max 



Vi - 

Ax ) 

v{x + Ax) — u(x) \ 1 

Tx — ) ix ' a ; 



v(x) — v(x — Ax) 

Ax 



(by definition (4.2)) 

< maxf^^ — cJ e fl n iti on (1.1) of p). 
xjiy\ x-y ) 


Thus (5.3) is proven. Let us now look at (5.4). We have 



Ax \ ( Ax 

~ ’*+') ~[ V ‘ + T 



(by definition (4.8)) 


Thus 

(5.5) 


6 i+ 1 = {vi+i-Vi ) — ^ max ^ 


AX f (Vi+2 - Vf+X V i+ i - Vi 


Ax 


Ax 


•HW. 4 4)) 


and therefore: 


c ^ , Ax(v i+l -Vi v i+l - 1\\ 

4 ‘+* * - ”<> - Tl Ax + — = 


0 


which enforces (5.4). This achieves the proof of (5.3-4) and therefore (4.10). So. 
the scheme is OSLC consistent. Let us now prove it is also second order accurate. 
To achieve that we shall prove: 

(5.6) f + 2 dx[E& t u{x) - E& t P&xv(x)] = O(Ax) 2 

for smooth functions u(x). 


We recall 


v{x) = Vi 4- 


(g ~ ft) 

Ax 


Si 


for x<_i < x < x t+ x, 
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where 



and S{ is defined by (4.8), (5.1). 

From the divergence theorem, this amounts to showing that: 

( 5 - 7 ) -£- t f Q [/(«(*<+*,*))- 

~~h j Q = 0{Ax) 2 

where u(x,t) and u&(x,t) are the solutions to (1.2) with initial data u(x) and v(x) 
respectively. 

This follows, if in each cell (x g -_x, x f+ x): 

(5.8a) u(x) — v(x) = c(x, x,)(Ax) 2 + O(Ax) 3 

with 

(5.8b) |c(x, x,) — c(x - Ax, x,_!)| < if Ax. 

But 

1 /**•'+ i 

u(x) - v(x) = ~r / [«(*) - u(y)\dy - (x - x,)s,- 

^ x 

= (x - x,)[u x (x.) - Si] + - ^(Ax) 2 ]« m (i.) + O(Ax ) 3 

Thus we must merely show: 

( 5 - 9 ) («*(*;) - «t) - (u x (x,--i) - s,_i) = 0((Ax) 2 ) 

but a simple calculation shows that 

tt*(zt) - Si = -^Ax\u xx (yi)\ +0(Ax) 2 


/(«*(*.•+*, ( 0 )]# 


( 5 . 10 ) 
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and we are finished. 

Remark. We do not recommend this convergent, fully second order scheme 
for practical calculations when shocks develop, because the limiter defined in (5.1) 
permits undershoots. This limiter agrees with the usual minmod limiter when a 
and b are both negative. Otherwise large jumps in the piecewise linear interpolant 
can be created and the differential equation solution operator may take some time 
to dissipate the resulting wiggles. See the calculations displayed in Section X below. 

For what we believe to be the state-of-the-art in practical non-oscillatory high 
order accurate shock capturing algorithms, see [6]. 

VI. A OSLC consistent method of lines and a second order order convergent Roe 
scheme. A MUSCL scheme as described above is not a practical scheme, becau e 
it involves the exact computation of (1.2), for 0 < t < At, when the initial value is 
a piecewise polynomial function. This task is, in general, impossible (except when 
the initial value is piecewise constant, as for the Godunov scheme). Actually, for 
schemes such as the ones described in Section IV (through definitions (4.7-8)) it is 
well known, that it is enough to know the exact solution only for x = z 1+ x and 
0 < t < At. This can be obtained through convenient approximations up to any- 
desired accuracy [1], [5]. 

A different approach consists of considering the semi-discrete scheme obtained 
when At — ► 0 (Az being fixed), the so-called method of lines. This leads to a 
stiff system of ordinary differential equations. Then one’s favorite ODE solver can 
be used to get a fully discrete scheme. The main advantage of this approach is 
the simplicity of the numerical expressions involved in the method of lines. As a 
matter of fact, the method of lines corresponding to any MUSCL scheme of the 
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form (4.7-8), can be written as follows (see [9]) 

( 6 . 1 ) 


diif 1 / . _ v 


( 6 . 2 ) 


f Ax Ax v 

fi+1 = h{Ui + —Si, U i+ 1 - —Si) 


where h(a,b) is the Godunov flux associated with /: 


(6.3) h(a, b) = /tGod(<b b) = e min^ef(a + s(b - a)) where e — sign(6 - a) 

In the special case of the OSLC consistent scheme described as above, we have 
(because of (4.8), (5.1)) 

(6.4) Si = ^ max(u, +1 - u,-, u,- - u,_i) 

Since the original scheme is OSLC consistent, its semi-discrete version is also OSLC 
consistent in the following sense: 

(6.5) tt. + ‘(0- MU < (p(u(0 0) -l + at) -l 

and therefore is fully entropy consistent. In the other hand, we have seen that the 
OSLC enforces: 


(6.7) 


Ax Ax 

l ” 3*+i n u% + 2 


(there is no increasing discontinuity in the piecewise linear approximations). 


Thus the Godunov flux hQ 0 d(a,b) is only used when the entries a, b satisfies 
a > b. But, we have the rather obvious result: 


LEMMA 2. Let f be a convex flux, ho 0 ^(a,b) the Godunov flux and hn oe (a,b) 
the Roe flux: 


(/(<■) •/</(») -/(«))(* -«> >0 
1/W •/(/(*) - /(«))(* - a) <0 


( 6 . 8 ) 
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Then for a > b, h Go ^{a, b) = h Roe (a, b ). 

Consequently, the OSLC consistent method of lines (6.1), (6.2), can be seen as 
the second order accurate MUSCL version of the Roe scheme, provided (6.4) is used 
to define the numerical slopes. This result is quite surprising since, as it has been 
often noted, the first order accurate Roe scheme is not (fully) entropy consistent. 

PROPOSITION 6. The semi-discrete second order MUSCL Roe scheme, de- 
fined by (6.4), is (strongly) OSLC consistent and therefore fully entropy consistent. 

VII. The Modified equation. It was pointed out to us by E. Harabetian, that 
the modified equation corresponding to our new scheme is of interest for the fol- 
lowing reason: On one hand, since the scheme is fully second order accurate, the 
perburbation term in the modified equation involves a third order derivative of the 
solution and therefore looks like a dispersive term (as in the KDV equation). On 
the other hand, since the OSLC is enforced, the perturbation term does not create 
wiggles (as in the KDV equation). After rather obvious computations, it turns out 
that: 


PROPOSITION 7. The modified equation associated with the second order Roe 
MUSCL scheme (6.1), (6.2), (6.4) is: 


(7.1) 

, U^ f(u) x "b (5|^xx|x — 0 

with 


(7.2) 

6 = 0( Ax 2 ) 


The study of (7.1), for <5 > 0 fixed, is in our opinion an interesting new open 
problem. Formally the solutions still satisfy the OSLC (1.9). But existence and 
uniqueness of solutions to (7.1) remain to be investigated. 
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They should be considered as limits of the smooth solutions of the fourth order 
parabolic equation: 

(7.2) Ut + f{u) x + <5|Wxx|x = £ u xxxxi £ 0 

A new entropy condition presumably follows from: 

(7.3) u xx has only increasing jumps. 

So we have the following: 

Conjecture. At least for smooth initial value, the Cauchy problem for (7.1) 
admits a unique solution u(t, x) such that: 

i) u xx (t,x) is in L{ oc 

ii) u xxx {t,x) > 7 > -oo, a.e. 

Moreover, u(t,x) satisfies the OSLC (1-9). 

This conjecture can be extended to: 

(7.4) u t +/(u)x - g{u xx ) x = 0 

for any concave function g (in (7.1), g(w) = — |ie|) 

Existence of travelling waves for (7.1), or (7.4), is a more feasible result, and is 
now under investigation. 

VIII. Discretization of the modified equation. The modified equation (7.1), and 
more generally equation (7.4), can be (semi-) discretized as follows. Let us define a 


15 


(first order) numerical flux g i+ i = g(wi, u>,+i) for the concave function g. It could 
be the Godunov flux: 


(8.1) 0God(M) = ■ g(a + s{b- a))], e= sgn(t-a), 

The E.O. flux (assuming g( 0) = <7'(0) = 0) 

(8.2) g EO {a,b) = g{b+)+g(a-) 

(where 6 + = max(o, 6), a_ = min(o, 6)) 


or the Lax Friedrichs flux of either form: 


(8.3) 


9hr(a, b) = g(- y^) - ^(6 - a) 

M + 9(b) X , L _a 
or ~(b~a) 

where the constant A satisfies : A > max|^(u)| 


Then, a method of lines approach can be designed for (7.4): 

(8.4) Ax ^t + A_ l/( ~ + Y ±l ') “ 6 9(^i, w.+i)] = 0 

(where A _a,- = a,- - a,_i) 


with: 

(8.5) Wi = (u,--i - 2 Ui + «,+ x) 


A fully discrete scheme can be designed by replacing by the Euler forward 
finite difference: 



Then we have, in the case g(w) = — |w|. 
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PROPOSITION 8. Scheme (8.6) satisfies the weak OSLC (2.3) under the sta- 
bility condition: 





+ ”?+! 


)l 


< 


< 


Ax 


Ax 2 6Af 


The proof is given in Appendix 2, in the case when the EO flux is used. 

As a consequence, when 6 = 0( Ax 2 ), (8.6) is a new OSLC consistent scheme 
for the conservation law (1.2). By a standard truncation analysis, this scheme is 
second order accurate in space (first order in time). Condition (8.7) is, then, a 
standard CFL condition. 


IX. Numerical results. Let us consider the Burgers equation: 


(9.1) 

with initial conditions: 

(9.2) 




if 0 < x < 4 
otherwise 


The solution at time t = 4 has been calculated exactly (Figure 1), by using 
Van Leer’s method (minmod slope limiter, Figure 2), or its “superbee” variant [11] 
of (Figure 3), and finally the OSLC consistent scheme described above (Figure 4). 
In each case, the time step was A t = 0.04 (25 time steps have been used), the mesh 
width Ax = 0.2. 

We see the excessive compression of the rarefaction wave using superbee and 
the expected undershoot at rarefactions using OSLC. 

Next we considered a linear equation 


(9,3) 


u t + u x = 0 
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with 


(9.4) « 0 (ar) = sin 3 a:, 

= 0 

We graph the solution as it advances one full period. 

Figure 5 shows the results using the minmod limiter, Figure 6 uses a compres- 
sive almost 3rd order TVD scheme devised in [14], and Figure 7 uses the OSLC 
scheme. The undershoot in this linear problem is more pronounced, as expected. 


7T. 


| z --!<0 


I w 

X ~ 2 >2 
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Appendix 1. (proof of Proposition 2) 


Let uo(x) be a piecewise smooth function defined for x 6 [0. L) and L-periodic. 
Let us denote by N the number of cells in such a way that: 


(All) 


L = N Ax 


u 0 has a finite number of discontinuity points £o < £i < £2 < • < £p- Let us 

introduce 


( A 1 .2) 


*, , (u 0 {x)-u 0 {y) 

p («o) = sup sup ( 

*=i p x V 


which is a finite constant (even if p(uq) = -Hoo). Assuming the propagation speed 
of the scheme to be uniformly bounded by C, let us define: 


( A 1 .3) 


7 n = {1 = 1, i V; dist(i'Ax, E) < CnAt + Ax} 


where E = {£ 0 , £ 1 , •••, £p}- Clearly, the values u” and do not depend on the 
values of u 0 around E if 1 £ I n . Therefore, using the OSLO consistency of the 
scheme, we get: 


( A 1 . 4) 


<±1 ~ < 


<p*{u 0 ) for i £ I n 


On the other hand, we always have; 


<±1 Z < < _1 


cmAf ’ 


since the scheme is OSLC consistent. 


Thus we have: 


(A1.6) 


u” +1 — u” < AxM”, where A/" — 


(anAt) 1 if i € I n 
p*{u 0 ) if i $ I n 
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let us now estimate the total variation at time nAt 


(A1.7) 
We have: 


TK = £ |«? +1 - <1 

TV n < £|u? +l - «" - AxM?\ + ^2 AxM? 


and, then (by (Al.6)): 


The second term of the RHS vanishes because of the periodic boundary conditions. 
Thus: 

(A1.8) 7V„ < 2 V AxM? = + 2 (AT - |/.|)A*p'(»o) 

by definition (Al.6). 


Clearly, by definition (Al.3), we have: 


(A1.9) 


i/»i< 


CnAt 

Ax 


+ 1 


Therefore: 

TV n <2-+2-^- + 2NAxp*(u 0 ) 
a anAt 

and thus: 


(Al.ll) TV n < 2- + 2-^- + 2Lp*(u 0 ) 

a a At 

since N Ax = L. 


Assuming Ax and At to be of the same order of magnitude, it follows that TV n 
is uniformly bounded. 

Thus, the scheme is total variation stable. Since it is also entropy consistent, 
the convergence the unique entropy solution follows through standard arguments. 
This completes the proof of Proposition 2. 
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Appendix 2 (Proof of Proposition 8) 


Let us consider scheme (8.6) in the case when 
(A2.1) g(w) = -|iu| = ~(w)+ - (-«->)+ (with a + = max(a,0)) 
and the E.O. flux is used: 


(A2.2) 


We get 


g(u,v) = 0 eo («, v ) = -v+ - (-«) + . 


(A2.3) 0 = - <) + (/(&&*) + + K-<) + ) 


1 


Let us introduce: 

(* !( ) rr +i = £«+,-<) 

we have (by definition (8.5)): 

1 




(A2.5) 

Using (A2.3), we get: 

A 7» 1 

(A2.6) 0 = — (P" + Y - P ( " } ) + Wi + «, where : 


(A2.7) 

and 


*i=A -(/( 


u ?e-l +"?^ \ _ ./■ «? + «&., 


M^)) 


(A2.8) 0* = «A.((< + ,) + + (-< +l ) + - (< +1 ) + - 

Let us rewrite 0a: 

(A2.9) ^ = «[(<+,) + + (-<+i) + - («? + i) + - (-<) + 

- (<+*)+ - (-<)++ (<) + + (-<->) J 


Thus 
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(A2.10) to > «[(-«&,) + - 2 (<m) + - 2 (~“") + + (<) J 

Let us now consider 0i : 

(A2.ll) = f{a) - 2/(6) 4- /(c) with 

(A2.12) a = + <+ 2 ), 6= \( u i +u i+i)> c =\( u ?-i +«*) 

Taylor expansions lead to: 

(A2.13) 0i =/(6) + f{b )(a - 6) + i?i - 2/(6) 

+ /(*)+ /'(*)(c - 6) + * 2 

where both f?i and i? 2 are positive, since / is convex. Therefore: 

(A2.14) 0i >/'(6)(a-26 + c) 

that is, using (A2.12) and definition (8.5) 

(A2.15) to > /'(i(<+< + ,))[< + < +1 ]^ 

Hence, (A2. 6-10-15) lead to: 

( A2.16) | * 7 (p?:i - P?+} ) < + <*»)*£) 

- £[(-*»)♦ -*(<•)+ - 2 H1 + («r) J 

Let us denote 


(A2.17) A = Pj+a, H = p, + i, C = Pi _i 

Using (A2.4), (A2.16) can be rewritten:. 


(A2.18) 
Ax 

a7 


(B" +1 - B*) < 


A” -C* 


• Ax 


Ax 2 


[(— .4 n + B) + - 2 (A n - B n ) + - 2{-B n + C n )+ + ( B n - C n )+ ] 
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That is: 

(A2.19) 

B n+ 1 < B n +v(A n -C n )-i[(-A n +B n ) + -2{A n -B n ) + -2{-B n +C n ) + +{B n -C n )+] 


where 


(A2.20) 


v — s 7 = a? 


Since: 


A n - C" = (A n - F n ) + (£ n - C n ) 


= (A n - F n )+ - (F n - A") + + (B n - C n )+ -(C 11 - B n )+ , 


we deduce: 


(A2.21) B n+l <B n + (A n - B n ) + (27 + v) 4- ( B n - A n ) + (-~i - v) 


and then: 


+ ( B n - C ")+(- 7 + */) + (<? n - B n ) + (27 - v) 


(A2.22) 


B n+i <B n 


Therefore if: 


+ [(/»*- B*)+ + (C* - B*)+] (27 + M) 
+ [(B"-^") + + (B"-C“) + ] (H-7) 


(A2.23) 
we get: 


M < 7 and 7 < - 


(A2.24) F n+1 < F" + ^[(A n - £ n ) + + ( C n - £")+] < max(A", B n , C n ) 
According to definitions (A2.17) and (A2.20), this means: 


(A2.25) 

provided: 

(A2.26) 


P* + \ l < max(^ + ,, /?_,) 


|W/ «"+«E 
Ax ' 2 


< —— < I 

Ax 2 Ax 6 


Since (A2.25) implies the weak OSLC condition, the proof of Proposition 8 is com- 
plete. 
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Fig. 2 MUSCL Scheme with 
minmod Limiter 
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Fig. 3 MUSCL Scheme with 
Superbee Limiter 
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LINEAR WAVE EQN . SOLUTION 
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Minmod Limiter 
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D 


Fig. 6 Compressive 3rd Order TVD Scheme 
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Fig. 7 OSLC Scheme 
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